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Abstract 

The production and decay of the X(3872) are analyzed under the assumption that the X is a 
weakly-bound molecule of the charm mesons D°D*° and D*°D°. The decays imply that the large 
D°D*° scattering length has an imaginary part. An effective field theory for particles with a large 
complex scattering length is used to derive factorization formulas for production rates and decay 
rates of X. If a partial width is calculated in a model with a particular value of the binding energy, 
the factorization formula can be used to extrapolate to other values of the binding energy and to 
take into account the width of the X. The factorization formulas relate the rates for production 
of X to those for production of D°D*° and D*°D° near threshold. They also imply that the line 
shape of X differs significantly from that of a Breit-Wigner resonance. 

PACS numbers: 12.38.-t, 12.39. St, 13.20. Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a narrow resonance near 3872 MeV discovered by the Belle collaboration 
in 2003 [1]. It has been observed through the exclusive decay B ± — > XK^ [1, 2] and through 
its inclusive production in proton- antiproton collisions [3, 4]. Its mass m x is extremely close 
to the threshold for the charm mesons D° and D*° [5]: 

m x - (m D o + m D *o) = +0.6 ± 1.1 MeV. (1) 

The upper bound on its decay width T x is [1] 

T x < 2.3 MeV (90% C.L.). (2) 

The X(3872) was discovered through its decay into J/ipn + n^ [1]. The Belle collaboration 
has recently presented evidence for the decays X — > J/ip it + it~it and X — > J/ip^ [6]. The 
decay X — > J/ipj implies that the X has positive charge conjugation. Upper bounds have 
been set on several other decay modes [7-11]. 

The nature of the X(3872) has not yet been definitively established. The presence of 
the J/ijj among its decay products motivates its interpretation as a charmonium state with 
constituents cc [12-14]. Two interpretations that are motivated by the proximity of m x to 
the D°D*° threshold are a hadronic molecule with constituents DD* [15-21] and a "cusp" 
associated with the D°D*° threshold [22, 23]. Other proposed interpretations include a 
tetraquark with constituents ccqq [24], a "hybrid charmonium" state with constituents ccg 
[25, 26], a glueball with constitutents ggg [27], and a diquark-antidiquark bound state with 
constituents cu + cu [28]. Measurements of the decays of the X can be used to determine 
its quantum numbers and narrow down the options [29-33]. The most predictive of the 
proposed interpretations are charmonium and DD* molecules. The C = — charmonium 
options are ruled out by the decay X — > J/tpj. Evidence ruling out or disfavoring each 
of the C = + charmonium options has been accumulating [5, 33, 34]. The most difficult 
charmonium state to rule out is the Xci(2P), partly because it can have resonant S-wave 
interactions with D°D*° and D*°D° that transform it into a DD* molecule [19]. 

The possibility that charm mesons might form molecular states was considered shortly 
after the discovery of charm [35-38] . The first quantitative study of the possibility of molec- 
ular states of charm mesons was carried out by Tornqvist in 1993 using a one-pion-exchange 
potential model. He found that the isospin-0 combinations of DD* and D*D could form 
weakly-bound states in the S-wave 1 ++ channel and in the P-wave 0~ + channel [39]. Since 
the binding energy is small compared to the 8.4 MeV splitting between the D°D*° thresh- 
old and the D + D*~ threshold, there are large isospin breaking effects [15, 16]. After the 
discovery of the X(3872), Swanson considered a potential model that includes both one- 
pion-exchange and quark exchange, and found that the C = + superposition of D°D*° and 
D*°D° could form a weakly-bound state in the S-wave 1 ++ channel [20]. Another mechanism 
for generating a DD* molecule is the accidental fine-tuning of the mass of the x cl (2P) to the 
D D*°/D*°D threshold which creates a DD* molecule with quantum numbers 1 ++ [19]. 

The assumption that the V(3872) is a weakly-bound DD* molecule is very predictive [19]. 
This assumption has been used by Voloshin to predict the rates and momentum distributions 
for the decays of X into D°D°7r° and D D ^ [17]. It has been used to calculate the rate for 
the exclusive decay of T(4S') into the X and two light hadrons [40], to estimate the decay 
rate for the discovery mode B + — > XK + [41], and to predict the suppression of the decay 
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rate for B° — > XK° [42]. The assumption that the X is a weakly-bound DD* molecule also 
has implications for inclusive production [43, 44]. 

In this paper, we point out that if X is a loosely-bound molecule, its short-distance decay 
rates and its exclusive production rates satisfy simple factorization formulas. In Section II, 
we illustrate the factorization formulas using a two-channel scattering model. In Section III, 
we show how the rate for a short-distance decay mode of the X, such as J/ip 7r + 7r~7r°, can 
be factorized into long-distance factor that involves the large scattering length and a short- 
distance factor that is insenstitive to a. In Section IV, we show how a production rate for 
the X, such as the decay rate for B + — > XK + , can be factorized into a long-distance factor 
and a short- distance factor. In Section V, we use factorization to calculate the shapes of 
the invariant mass distributions of D°D*° and D*°D° near threshold. In Section VI, we use 
factorization to calculate the line shape of the X in any of its short- distance decay modes. 
The line shape depends on the real and imaginary parts of the scattering length and can 
differ substantially from a conventional Breit-Wigner resonance. A summary of our results 
is given in Section VII. 



II. FACTORIZATION IN A SIMPLE SCATTERING MODEL 
A. Universality for large scattering length 

Nonrelativistic few-body systems with short-range interactions and a large scattering 
length a have universal properties that depend on the scattering length but are otherwise 
insensitive to details at distances small compared to a [45]. In any specific system, there is a 
natural momentum scale A that sets the scale of most low-energy scattering parameters. The 
scattering length is large if it satisfies \a\ ^> A -1 . Universality predicts that the T-matrix 
element for 2-body elastic scattering with relative momentum p <C A is 

T(p) = -MH^ t (3) 
— 1/a — ip 

where /i is the reduced mass of the two particles. If a is real and positive, universality 
predicts that there is a weakly-bound state with binding energy 

E * = < 4) 

The universal momentum-space wavefunction of this bound state is 

= 2,1/2 - 5 

p l + l/cr 

The universal amplitude for transitions from the bound state to a scattering state consisting 
of two particles with relative momentum p <C A is 



A x = ^a-W (6) 

These results are all encoded in the universal expression for the truncated connected tran- 
sition amplitude: 

a ( e) = - ( 7) 
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FIG. 1: The geometric series of Feynman diagrams whose sum is the universal amplitude A(E). 

The universal amplitude in Eq. (7) can be obtained from a local effective field theory 
for the two particles. The particles interact through an S-wave contact interaction with 
Feynman rule — i(2n/ fi)a (A), where the parameter a (A) is a bare scattering length that 
depends on the ultraviolet momentum cutoff A. The amplitude A(E) can be obtained by 
summing the geometric series of Feynman diagrams in Fig. 1: 

A(F] _ (27r//x)a (A) 
1 ] l + (27r//z)a (A)L(A,£)' 1 ] 

where L(A, E) is the amplitude for the propagation of the particles between successive 
contact interactions: 

L(A,£0 = ^(A-|v C 2^)- (9) 

Renormalization is accomplished by eliminating ao(A) in favor of the physical scattering 
length: 

ao(A) (10) 



l + (2/7r)Ao (A)' 

With this substitution, the expression for A(E) in Eq. (8) reduces without approximation 
to the universal result in Eq. (7). Note that the scattering length a can be tuned to ±oo by 
tuning the bare scattering length to a critical value of order A -1 : 

oo(A) — -\A-\ (11) 

If the 2-body system has inelastic scattering channels, the large scattering length a will 
have a negative imaginary part. It is convenient to express the complex scattering length in 
the form 

- = 7re + r/im, (12) 
(X 

where 7 re and 7; m are real and 7i m > 0. In the case 7 re > where there is a weakly-bound 
state, it can decay into the inelastic channel. The expression for the binding energy on the 
right side of Eq. (4) is complex- valued. Its real part E x , po ie and its imaginary part T^/2 are 
given by 

£x, P oic = (7r C -7im)/(2/i), (13a) 
= 27 rc 7 im //i. (13b) 
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These quantities specify that there is a pole in the S-matrix at the energy E — — .Ex, pole — 
iTx/2. As we shall see in Section VI, Tx can be interpreted as the full width at half- 
maximum of a resonance in the inelastic channel provided 7; m < 7 rc . The peak of the 
resonance is below the threshold for the two particles by 

E x = 7 r V(2/i). (14) 
We therefore interpret E x as the binding energy rather than Ex,po\c- 



B. Two- channel model 

Cohen, Gelman, and van Kolck have constructed a renormalizable effective field theory 
that describes two scattering channels with S-wave contact interactions [46]. We will refer 
to this model as the two- channel scattering model. An essentially equivalent model has been 
used to describe the effects of A A states on the two-nucleon system [47]. The parameters 
of this model can be tuned to produce a large scattering length in the lower energy channel. 
It can therefore be used as a simple model for the effects on the D°D* /D*°D system of 
other hadronic channels with nearby thresholds, such as J/ipp, J/ipcu, and D ± D* =F . 

The two-channel model of Ref. [46] describes two scattering channels with S-wave contact 
interactions only. We label the particles in the first channel la and lb and those in the 
second channel 2a and 2b. We denote the reduced masses in the two channels by /ii and 
fi2- Renormalized observables in the 2-body sector are expressed in terms of 4 parameters: 
three interaction parameters an, ai 2 , and a 22 with dimensions of length and the energy gap 
A between the two scattering channels, which is determined by the masses of the particles: 

A = m 2a + m 2b - (m la + m lb ). (15) 

The scattering parameters in Ref. [46] were defined in such a way that an and a 22 reduce in 
the limit ai 2 — > ±oo to the scattering lengths for the two channels. The truncated connected 
transition amplitude A(E) for this coupled-channel system is a 2 x 2 matrix that depends 
on the energy E in the center-of-mass frame. If that energy is measured relative to the 
threshold m\ a + m\ b for the first scattering channel, the inverse of the matrix A(E) is 1 

" 2tt V y/mti/a 12 /i 2 [- l/a 22 + v/2/i 2 (A-E)] 

The square roots are defined for negative real arguments by the prescription E — > E + ie 
with e — > + . The amplitudes defined by Eq. (16) are for transitions between states with the 
standard nonrelativistic normalizations. The transitions between states with the standard 
relativistic normalizations are obtained by multiplying by a factor y / 2m^ for every particle 
in the initial and final state. We will need explicit expression for the 11 and 12 entries of 



1 The expression for the matrix T s 1 in Eq. (2.18) of Ref. [46] should be equal to A(E) 1 evaluated at 
E = p 2 /(2m). There is an error in the 22 component of T s 1 : the square root \J p 2 — 2/i2A should be 
V / (M2//^i)p 2 - 2^ 2 A. 



(16) 
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this matrix: 



2tt 



An(E) = — + - -= l/a 22 + y/2fi 2 (A - E) 



A 12 (E) 



Hi \ an 



2tt 



12 



-1 



- a 12 



— - + y/-2\i x E — — + a/2/x 2 (A - E) J 
an J [ a 22 J/ 



-i 



(17a) 



.(17b) 



The T-matrix element for the elastic scattering of particles in the first channel with 
relative momentum p is obtained by evaluating Au(E) in Eq. (17a) at the energy E = 
P 2 /(2/ii): 



Tn(p) = 



2tt 



ip - — 

ai 



~ 1/022 + A/2yU 2 A - (/i 2 //il)p 2 



-1 



A*i \ fl n "i 2 
Setting T n (0) = — 27r/(/i 1 a), we can read off the inverse scattering length 1/a: 



(18) 



a an af 2 



a/2/i 2 A - l/a 22 



(19) 



If there is a bound state with energy — /t 2 /(2/ii) below the scattering threshold for the 
first channel, the matrix A(E) given by Eq. (16) has a pole at E — —K 2 /(2fii). The binding 
momentum k satisfies 



K— h — 

an a 12 



l/a 22 + a/2u 2 A + (/i 2 //ii)/t 2 



(20) 



The behavior of the matrix ^4(-E) as the energy E 1 approaches the pole associated with the 
bound state is 



A(E) 



1 



E+K 2 /(2fi 



(21) 



The components Axi and ^4^2 of the column vector are the amplitudes for transitions from 
the bound state to particles in the first and second channels, respectively. The column vector 
is an eigenvector of the matrix A(E)^ 1 in Eq. (16) with eigenvalue zero, so its components 
must satisfy 



ui[-l/a n + re] ^4xi + [y/fJiitJv/avi] A X 2 = 0. 



(22) 



C. Two-channel model with large scattering length 

The two-channel model of Ref. [46] can be used as a phenomenological model for a system 
with a large scattering length a in the first channel. The large scattering length requires 
a fine-tuning of the parameters an, a 22 , ai 2 , and A. There are various ways to tune the 
parameters so that a — > ±oo. For example, if an < a 2 2 /a 22 , the energy gap A can be tuned 
to the critical value where the right side of Eq. (19) vanishes. Alternatively, the scattering 
parameter an can be tuned to the critical value — a 2 2 [y / 27t 2 A~ — l/a 22 ]. The coefficients in 
the low-energy expansion of 7ii(p) _1 are proportional to various powers of 1/an, l/ai 2 , and 
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s/2jI^K. We assume that these momentum scales are comparable in magnitude. We refer 
to that common momentum scale as the natural low-energy scale and we denote it by A. 

For \a\ 3> A -1 and \E\ <C A 2 /(2/ii), the amplitude Au(E) in Eq. (17a) approaches the 
universal amplitude A(E) in Eq. (7) with ji = It follows that for p <C A the T-matrix 
element %_i(p) in Eq. (18) approaches the universal T-matrix element T(p) in Eq. (3). 
For a ^> A -1 , the solution to Eq. (20) for the binding momentum k approaches 1/a, so 
/t 2 /(2/i) approaches the universal binding energy E x in Eq. (4). Finally the amplitude Axi 
for transitions from the bound state to particles in the first channel, which is defined in 
Eq. (21), approaches the universal amplitude Ax in Eq. (6). 

There are also universal features associated with transitions to the second channel. If 
\a\ ^> A -1 and \E\ <C A 2 //ii, the leading term in the transition amplitude Au(E) in 
Eq. (17b) reduces to 

Al2 (E) = -^3j*l [y / 2^A - l/a 22 1 ~\{E), (23) 
ai2 L J 

where A(E) is the universal amplitude in Eq. (7) with [i replaced by For a ^> A^ 1 , the 
leading term in the amplitude Ax2 for transitions of the weakly-bound state X to particles 
in the second channel, which is defined in Eq. (21), reduces to 

V'/'I '/' • — - 

•A-X2 — 



Ol2 



V / 2^A-l/a 22 A x , (24) 



where Ax is the universal amplitude in Eq. (6). Note that the ratio Au(E)/Ax2 of the 
amplitudes in Eqs. (23) and (24) is a universal function of a and E only. 

The expressions for An{E) and Ax2 in Eqs. (23) and (24) are examples of factorization 
formulas. They express the leading terms in the amplitudes as products of the same short- 
distance factor and different long-distance factors A(E) and Ax- The long-distance factors 
involve the large scattering length a. The limit \a\ — > oo has been taken in the short- 
distance factors. The conditions \E\ <C A 2 /(2/i) or E = —Ex require the particles in the 
second channel to be off the energy shell by approximately A. In the short time 1/A allowed 
by the uncertainty principle, those particles can propagate only over short distances of order 
(2/x 2 A) -1 / 2 . This is small compared to the distance scales (2n\E\Y l l 2 or \a\ associated with 
the particles in the first channel. Thus as far as they are concerned, the particles in the 
second channel act only as a point source for particles in the first channel. The amplitudes for 
particles from such a point source to evolve into particles of energy E and into the weakly- 
bound state are L(A, E)A(E) and L(A, — E x )Ax, respectively. By using the conditions 
\E\,E X <C A 2 /(2/x), these amplitudes reduce to (fiA/-K 2 )A(E) and (fiA/7T 2 )Ax, respectively. 
In these expressions, the short-distance factors are identical and the long-distance factors 
are the same as those in Eqs. (23) and (24). 



D. Unstable particle in the second channel 

Now let us suppose one of the scattering particles in the second channel has a nonzero 
width. We take that particle to be 2b. We assume that its width T 2 b arises from its decay 
into particles with relativistic momenta that are much greater than the ultraviolet cutoff 
Auv that defines the domain of validity of the two-channel model. The momenta of the 
decay products are therefore also much greater than a/2/z 2 A. We assume that T 2 b is small 



7 



compared to the mass m 2 b, but not necessarily small compared to the energy gap A between 
the two channels. This makes it necessary to take into account the contribution to the 
self-energy of particle 2b from the coupling to its decay products. 

Taking into account the self-energy of particle 2b would modify the term ^2{i 2 (A — E) 
in the inverse of the matrix of transition amplitudes given in Eq. (16). That term arises 
from the amplitude for the propagation of particles in the second channel between contact 
interactions, which is given by the integral 

/ Wf E-A-^),,, = % ( A - - i^MA . (25) 

The cutoff constrains the momentum to the region in which the nonrelativistic approximation 
for the energy of the particle 26 is valid. In this region, the self-energy II can be expressed 
as a function of E' = E — A — p 2 /(2ji 2 ). It can be taken into account by replacing A in 
the integral in Eq. (25) by A + H(E'). The assumption that the decay products of particle 
2b have relativistic momenta comparable to m 2b implies that their contributions to H(E') 
have significant dependence on E' only for variations in E' that are comparable to m 2 b. 
For energies satisfying \E\ <C Auv/(2/i2) and loop momenta p < Auv, the dependence on 
E' can be neglected and the argument of n(£") can be set to a constant, such as —A. 
The prescription for taking into account the self-energy then reduces to replacing A in the 
integral in Eq. (25) by A + Il(— A). The real part of Il(— A) can be absorbed into A so that 
it becomes the physical threshold. The imaginary part of II(— A) is related to the width of 
particle 2b: Imll(— A) = — T 2 b/2. Thus the leading effect of the self-energy can be taken 
into account by replacing A in Eq. (16) by the complex-valued energy gap 

A = m 2a + m 2b - (m la + m lb ) - iT 2b /2. (26) 

If A is complex, the solution to Eq. (20) for the binding momentum k is complex. It 
determines the pole mass mx, po ic and the width Tx of the weakly-bound state according to 

m la + m u - k 2 / (2/ii) = m x , P oic - i^x/2. (27) 

The imaginary part reflects the fact that the bound state can decay into particle 2a and 
decay products of particle 2b. The quantities m x , vo \c and Y x in Eq. (27) give the location 
of a pole in the S-matrix. They need not have the standard interpretations as the location 
of the peak and the full width at half maximum of a Breit-Wigner resonance. 



E. The D D*°/D*°D System 

The energy difference between the mass of X and the D°D*° threshold, which is given 
in Eq. (1), is small compared to the natural energy scale for binding by the pion exchange 
interaction: m 2 /2/j, ~ 10 MeV, where \x is the reduced mass of D° and D*°. The unnaturally 
small value of the energy difference implies that if the X couples to D°D*° and D*°D°, the 
S-wave scattering lengths for those channels must be large compared to the natural length 
scale 1/rrin associated with the pion exchange interaction. We assume that there is a large 
scattering length a in the channel (DD*)° + with even charge conjugation defined by 

\(DD*)° + ) = -i= (\D°D*°) + \D*°D )) . (28) 
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If the scattering length in the C = — channel is negligible in comparison, the scattering 
lengths for elastic D°D*° scattering and elastic D*°D° scattering are both a/2. We identify 
the X as a bound state in the (DD*)° + channel. 

The decays of the X imply that the scattering length a is complex-valued. It can be 
parameterized in terms of the real and imaginary parts of 1/a as in Eq. (12). Our inter- 
pretation of X as a bound state requires 7 re > 0. The energy difference in Eq. (1) puts an 
upper bound on 7 re : 

7re < 40 MeV (90% C.L.). (29) 
The upper bound on the width in Eq. (2) puts an upper bound on the product of 7 re and 

Tim- 

7rc7im < (33 MeV) 2 (90% C.L.). (30) 

There is also a lower bound on the width of the X from its decays into D°D°it and D°D°'y, 
which both proceed through the decay of a constituent D*. These decays involve interesting 
interference effects, but the decay rates have smooth limits as the binding energy is tuned 
to [17]. In this limit, the constructive interference increases the decay rate by a factor of 
2: the partial width of X reduces to 2 r[L>* ]. The width of D*° has not been measured, but 
it can be deduced from other information about the decays of D*° and D* + . Using the total 
width of the D* + , its branching fraction into D + tt°, and isospin symmetry, we can deduce 
the partial width of D*° into D°n°: T[D*° -> D tt ] = 42 ± 10 keV. The total width of the 
D*° can then be obtained by dividing by its branching fraction into D°ir°: T[D*°] = 68 ± 16 
keV. The sum of the partial widths of X into D°D n° and D°D°-f is therefore 136 ± 32 keV. 
The resulting lower bound on the product of 7 re and 7 im is 

7 rc7im > (7 MeV) 2 (90% C.L.). (31) 

By combining this with the upper bound on 7 re in Eq. (29), we can infer that 7 im > 1 MeV. 

III. SHORT-DISTANCE DECAYS OF X 

The decay modes of the X(3872) can be classified into long-distance decays and short- 
distance decays. The long-distance decay modes are D°D°ir and _D°D°7, which proceed 
through the decay of a constituent D*° or D*°. These decays are dominated by a component 
of the wavefunction of the X in which the separation of the D and D* is of order l/|a|. These 
long-distance decays involve interesting interference effects between the D°D*° and D*°D° 
components of the wavefunction [17]. The short- distance decays involve a component of the 
wavefunction in which the separation of the D and D* is of order l/m n or smaller. Examples 
are the observed decay modes J ftp 7r + 7r~, J/ip 7r + 7r~7r°, and J/ip'y. 

Short-distance decays of the X into a hadronic final state H involve well-separated mo- 
mentum scales. The DD* wavefunction of the X involves the momentum scale l/|a| set 
by the large scattering length. The transition of the DD* to H involves momentum scales 
m w or larger. We will refer to momentum scales of order l/|a| and smaller as long-distance 
scales and momentum scales of order m w and larger as short- distance scales. We denote the 
arbitrary boundary between these two momentum regions by A. 
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The separation of scales \a\ ^> l/m w in the decay process X — > H can be exploited 
through a factorization formula for the T-matrix element: 

T[X - if] = A x x A shoit [(DD*)° + - H}. (32) 

In the long-distance factor, Ax is the universal amplitude given in Eq. (6) and the factor of 
\J1mx takes into account the difference between the standard nonrelativistic and relativistic 
normalizations of states. If the complex scattering length is parameterized as in Eq. (12), 
this factor is 



A x = — ^ (7rc + ilim) 1/2 • (33) 

The short-distance factor *4. s hort in Eq. (32) is insensitive to a, and one can therefore take the 
limit \a\ — > oo in this factor. The factorization formula in Eq. (32) can serve as a definition 
of the short-distance factor. The content of the factorization statement then resides in the 
fact that, up to corrections suppressed by powers of l/(am n ), the same short- distance factor 
appears in the factorization formula for the T-matrix element for the scattering process 
D°D*° — > if at energies E near the D°D*° threshold: 

T[D°D*° ^ H} = ^-= V4m D om D .o A(E) x A short [(DD*)° + -> H\. (34) 
v2 

In the long-distance factor, the \/y/2 is the amplitude for D°D*° to be in the channel (DD*)° + 
with the large scattering length, the factor of y^moom^o takes into account the difference 
between the standard nonrelativistic and relativistic normalizations of states, and A(E) is 
the universal amplitude given in Eq. (7). If the complex scattering length is parameterized 
as in Eq. (12), this factor is 

A(E) = piL (35) 

~7rc - H7im + \/2\lE) 

The factorization formulas in Eqs. (34) and (32) are analogous to those in Eqs. (23) and 
(24) for the two-channel model with a large scattering length in the first channel. 

The factorization formulas in Eqs. (32) and (34) can be motivated diagrammatically by 
separating virtual particles into soft particles and hard particles according to whether they 
are off their energy shells by less than or by more than A 2 /(2/i), where A is the arbitrary 
momentum separating the long-distance scale l/\a\ and the short- distance scale m n . Any 
contribution from soft particles inside a subdiagram all of whose external legs are hard can 
be Taylor-expanded in the momentum of the soft particles, leading to suppression factors of 
l/(aA). The diagrams with the fewest suppression factors will be ones that can be separated 
into a part for which all the internal lines are hard particles and a part that involves only 
soft particles. This separation leads to the factorization formula. 

The leading terms in the T-matrix element for X — > H in the limit |a|m w ^> 1 can be 
represented by the Feynman diagrams in Fig. 2 and can be expressed as 

T[X ^H] = V2^x~ J ^- 3 HP)A^[(DD*)1 - H], (36) 

where ip(p) is the universal wavefunction in Eq. (5). The factor A^, which is represented 
by a dot in Fig. 2, is an amplitude for the transition (DD*)° + — > H in which all virtual 
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FIG. 2: Feynman diagram for X — ► J/-0 7r + 7r that scale like (am n )~ - 1 / 2 . The (DD*)° + wavefunction 
of the X is the integral over the loop energy of the product of the blob and the attached propagators. 




FIG. 3: Feynman diagrams for D°D*° — ► J/i[)it + it that are enhanced near the D°D*° threshold 
by a factor of am n . The blob represents the geometric series of Feynman diagrams in Fig. 1. 

particles are off their energy shells by more than A 2 /(2/i). It is therefore insensitive to the 
relative momentum p of the D and D*. If that momentum dependence is neglected and 
if the integral in Eq. (36) is regularized by a momentum cutoff \p\ < A, the wavefunction 
factor reduces in the limit \a\ ^> A -1 to 




The factorization formula in Eq. (32) is then obtained by absorbing a factor of (///7r 2 )A into 
A^ to obtain the short-distance factor: 

A shoit [(DD*)° + - H\ = (^A) AW[{DD% - H\. (38) 

Since the T-matrix element in Eq. (32) is independent of the arbitrary separation scale, the 
dependence on A must cancel between the two factors on the right side of Eq. (38). 

The leading term in the T-matrix element for D°D*° — > H in the limits \a\ ^> A -1 and 
E -C A 2 / (2/z) can be represented by the Feynman diagram in Fig. 3 and can be expressed 
as 

T[D°D*° ^H] = ^= y /4m DO m D .o A{E) L(A, E) A {K) [{DD% -> H]. (39) 
v2 

The factor L(A, E) is the amplitude for the propagation of the D and D* between successive 
contact interactions, which is given in Eq. (9). The approximation E A 2 /(2/i) justifies 
neglecting the \f2jlE term in L(A, E). The factorization formula in Eq. (34) is then obtained 
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by absorbing the remaining term (/x/7r 2 )A into „4.( A ) to obtain the short-distance factor in 



The factorization formula for the T-matrix element in Eq. (32) implies a factorization 
formula for the decay rate for X — > H. The decay rate T[X — > H] can be expressed as the 
product of a short-distance factor and the long-distance factor 



Using the expressions for the binding energy and the total width of the X in Eqs. (14) and 
(13b), the long-distance factor in Eq. (40) can be expressed as 



If the partial width for a short-distance decay mode of the X has been calculated using 
a model with a specific binding energy, the factorization formula for the decay rate can be 
used to extrapolate the prediction to other values of the binding energy and to take into 
account the effect of the width of the X. This is useful because numerical calculations in 
models often become increasingly unstable as the binding energy is tuned to zero. Swanson 
has estimated the partial widths for various short- distance decays of X using a potential 
model, but only for binding energies down to about 1 MeV and without taking into account 
the effect of the width of the X [20, 21]. His predictions can be extrapolated to other 
values of the binding energy and the width of the X can be taken into account by using the 
long-distance factor in Eq. (41). 

IV. PRODUCTION OF X 

The production of X necessarily involves the long-distance momentum scale l/|a| through 
the (DD*)° + wavefunction of the X. The production also involves much larger momentum 
scales. Unless there are already hadrons in the initial state containing a c and c, the produc- 
tion process involves the scale m c associated with the creation of a cc pair. Even if the initial 
state includes hadrons that contain c and c, such as J/ip or D + and D~ , the production 
process involves the scale m n associated with the formation of the D° and D*° that bind 
to form the X. We will define a short- distance production process to be one for which the 
initial state either does not include any of the charm mesons D°, D*°, D*°, or D°, or if 
it does, the momentum of the charm meson in the rest frame of the X is of order or 
larger. All practical production mechanisms for X in high energy physics are short- distance 
processes. Long-distance production mechanisms could arise in a hadronic medium that 
includes charm mesons, such as that produced by relativistic heavy-ion collisions. 

In a short-distance production process, the separation between the long-distance scale 
l/\a\ and all the shorter-distance momentum scales can be exploited through a factorization 
formula that expresses the leading term in the production rate as the product of a long- 
distance factor that involves a and a short- distance factor that is insensitive to a. To be 
definite, we will consider the specific production process B — > XK. The factorization for 
any other short-distance production process will have the same long-distance factor but a 
different short-distance factor. 



Eq. (38). 




(40) 




(41) 
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FIG. 4: Feynman diagram for B — > XK that scales like (am n ) 1 ' 2 . The (DD*) ^ wavefunction of 
the X is the integral over the loop energy of the product of the blob and the attached propagators. 

There are many momentum scales that play an important role in the decay B — > XK, 
ranging from the extremely short-distance scales mw and nib associated with the quark decay 
process b — > ccs to the smaller short- distance scales Aqcd and rn^ involved in formation of the 
final-state hadrons to the long-distance scale l/\a\ associated with the (DD*)° + wavefunction 
of the X. We denote the arbitrary boundary between the long-distance momentum region 
and the short- distance momentum region by A. 

The separation between the long-distance scale l/\a\ and all the short-distance momen- 
tum scales in the decay B — > XK can be exploited through a factorization formula for the 
T- matrix element: 



In the long-distance factor, Ax is the universal amplitude in Eq. (33). The short- distance 
factor in Eq. (42) is insensitive to a and one can therefore take the limit \a\ — > oo in that 
factor. The factorization formula in (42) can serve as the definition of the short- distance 
factor. The content of the factorization statement then resides in the fact that the same 
short-distance factor appears in the factorization formula for the T-matrix element for the 
decay B — > D°D*°K when the DD* invariant mass is near the D°D*° threshold. This 
factorization formula is discussed in Section V. 

The factorization formula in Eq. (42) can be motivated diagrammatically by separating 
the loop integrals in the decay amplitude according to whether the virtual particles are off 
their energy shells by less than or by more than A 2 /(2/i). The leading terms in the decay 
amplitude for B — > XK are suppressed only by a factor of (am,)" 1 / 2 . These terms can be 
represented by the Feynman diagram in Fig. 4 and can be expressed in the form 



where ip(p) is the universal wavefunction in Eq. (5). The factor A A **, which is represented 
by a dot in Fig. 4, is an amplitude for the decay B — > (DD*)° + K in which all virtual particles 
are off their energy shells by more than A 2 /(2/i). It is therefore insensitive to the relative 
momentum p of the D and D*. If that momentum dependence is neglected, the wavefunction 



T[B -> XK] = Abort [5 -> (DD*)° + K] x AxV^- 



(42) 




(43) 
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factor in Eq. (43) reduces to Eq. (37). The factorization formula in Eq. (42) then requires 
the short- distance factor to be 

A short [B -> {DD*f + K\ = AW[B -> {DD*f + K} (£ A ) . (44) 

Since the T-matrix element in Eq. (42) is independent of the arbitrary separation scale, the 
dependence on A must cancel between the two factors on the right side of Eq. (44). 

We proceed to use the factored expression in Eq. (42) to evaluate the decay rate for 
B + — > XK + . Lorentz invariance constrains the short-distance amplitude «4 s hort at the DD* 
threshold to have the form 

Abort [B + -> (DD*)° + K + ] = c + P-e D *, (45) 

where P is the 4-momentum of the B + and e^* is the polarization 4- vector of the D*. Heavy 
quark spin symmetry guarantees that the polarization vector to* can be identified with the 
polarization vector ex of the X. The decay rate is obtained by squaring the amplitude in 
Eq. (42), summing over the spin of the X, and integrating over phase space. The resulting 
expression for the decay rate is 

r[B+ - XK + ] = \c + \* XV2{mB ^ mK) \A x \\ (46) 

where X(x,y,z) is the triangle function: 

\{x, y, z)=x i + y 4 + z i - 2{x 2 y 2 + y 2 z 2 + z 2 x 2 ). (47) 

The long-distance factor \Ax\ 2 is given in Eq. (40). The result in Eq. (46) was obtained in 
Ref. [41] for the special case 7 im = and used to estimate the order of magnitude of the 
decay rate for B + — > XK + . The estimate is consistent with the measurement of the product 
of the branching fractions for B + — > XK + and X — > <//■(/> 7r + 7r~ [1] provided J/ifjir + TT^ is 
one of the major decay modes of X. 

The factorization formula for the decay rate for B° — > XK° has the same form as in 
Eq. (46) except that the coefficient c + in the short-distance decay amplitude in Eq. (45) 
has a different value. In Ref. [42], it was pointed out that the decay rate for B° — > XK° 
should be suppressed compared to B + — > XK + . That suppression can be understood by 
considering the short-distance amplitude for B — > DD*K. The dominant contributions to 
most decay amplitudes of the B meson are believed to be factorizable into the product 
of matrix elements of currents. The factorizable contributions to the decay amplitude for 
B + — > (DD*)°^K + have three terms: the product of B + — > D*° and — > D°K + matrix 
elements, where is the QCD vacuum, the product of B + — > D° and — > D*°K + matrix 
elements, and the product of B + — > K + and — > (DD*)° + matrix elements. The factorizable 
contributions to the decay amplitude for B° — > (DD*)°^K have only one term: the product 
of B° —> K° and — > (DD*) ^ matrix elements. Heavy quark symmetry implies that the 
-> (DD*)° + _ matrix element vanishes at the D°D*° threshold. The decay B° -> (DD*)° + K° 
near the D°D*° threshold must therefore proceed through nonfactorizable terms in the decay 
amplitude. The resulting suppression of the coefficient c + in the short- distance factor for 
B° — > (DD*)° + K° results in a suppression of the rate for B° — > XK° relative to the rate for 
B + — > XK + . In Ref. [42], a quantitative analysis of Babar data on the branching fractions 
for B -> D^D^K [48] was used to estimate the suppression factor to be an order of 
magnitude or more. 
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D° 




FIG. 5: Feynman diagrams for B — > DDK that are enhanced near the D D threshold by a 
factor of am-n- The blob represents the geometric series of diagrams shown in Fig. 1. 

V. PRODUCTION OF D°D*° NEAR THRESHOLD 

It was pointed out in Ref. [41] that the identification of X as a DD* molecule could be 
confirmed by observing a peak in the invariant mass distribution for D°D*° (or D*°D°) near 
the DD* threshold in the decay B — > DD*K. The shape of that invariant mass distribution 
was given for a real scattering length a. The shape would be the same for any other short- 
distance production process. In this section, we consider the effect of an imaginary part of 
the scattering length on the DD* invariant mass distribution for a short-distance production 
process. To be specific, we consider the short-distance production process B — > D°D*°K. 

The separation between the long-distance scale l/|a| and all the short-distance momen- 
tum scales in the decay B — > D°D*°K can be exploited through a factorization formula for 
the T-matrix element: 

T[B -> D°D*°K] = A shoit [B -> (DD*)° + K] x A{E) V4m D om D .o — . (48) 

v2 

In the long-distance factor, A{E) is the universal amplitude in Eq. (35) and the factor l/\/2 
is the amplitude for D°D*° to be in the channel (DD*)° + with the large scattering length. 
The short- distance factor ^4 s h or t is the same as in the factorization formula for B — > XK in 
Eq. (42). 

The factorization formula in Eq. (48) can be motivated diagrammatically by separating 
the loop integrals in the decay amplitude according to whether the virtual particles are off 
their energy shells by less than or by more than A 2 /(2yu). There are terms in the T-matrix 
element for the decay B — > D°D*°K that are enhanced near the DD* threshold by a factor 
of am-n- These terms can be represented by the Feynman diagrams in Fig. 5 and can be 
expressed in the form 

T[B -> D°D*°K} = A w [B -> {DD*f + K} L(A, E) A(E) V4m D om D ,o —=. (49) 

V2 

The first factor A^ A \ which is represented by the dot in Fig. 5, is an amplitude for the 
decay into (DD*)°^K in which all virtual particles are off their energy shells by more than 
A 2 /(2/i). It is therefore insensitive to the relative momentum p. The second factor L(A, E) 
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FIG. 6: The DD* invariant mass distribution in B — > D°D*°K for y im = 10 MeV and various 
values of |7re|/7im- The horizontal axis is the difference E = M — (mjjo + m D *o) between the 
invariant mass M and the D°D*° threshold. 



is the amplitude for the propagation of the D and D* between contact interactions, which 
is given in Eq. (9). The condition \E\ <C A 2 /(2/i) implies that the ^—2fiE term in L(A,E) 
can be neglected. The resulting expression for the T-matrix element is the factorization 
formula in Eq. (48), with the short- distance factor A^ort given in Eq. (44). 

We proceed to use the factorized expression for the decay amplitude in Eq. (48) to 
calculate the DD* invariant mass distribution near the D°D*° threshold in the decay B + — > 
D°D*°K + . Lorentz invariance constrains the short- distance amplitude „4 s hort at the DD* 
threshold to have the form in Eq. (45). The decay rate is obtained by squaring the amplitude 
in Eq. (48), summing over the spins of the D*°, and integrating over phase space. The 
resulting expression for the differential decay rate with respect to the DD* invariant mass 
M near the D°D*° threshold is 



dT 

dM 



[B + -> D°D*°K + ] = |c+| 



,jiX 3/2 (m B ,M, m K ) 
2567r 3 m|M 2 



A X / 2 (M, 



m D o,m D *o)\A(E)\ 2 , (50) 



where E is the energy of the D°D*° in its rest frame relative to the D°D*° threshold: 

E = M - (m D o +m D *o). (51) 

We have used the fact that M is near the D°D*° threshold to replace a factor of m D om D *o in 
Eq. (50) by fiM. The result in Eq. (50) was obtained previously in Ref. [41] for the special 
case 7i m = 0. For M near the D°D*° threshold, the only significant variation with M is 
through the long-distance factor |^l(i?)| 2 and the threshold factor 



) w 1M\[2ixE. (52) 
If the complex scattering length is parameterized as in Eq. (12), the long-distance factor is 

47T 2 //i 2 



\A{E)Y 



((2^)V2 +7im )2 + 7] 



(53) 
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The shape of the D°D*° invariant mass distribution in Eq. (50) is given by the factor 
y/2jlE \ A(E)\ 2 . Note that it depends on 7 rc and 7i m but not on the sign of 7 rc . The invariant 
mass distribution is shown in Fig. 6 for 7 im = 10 MeV and three values of |7 re |: 10, 20, 
and 30 MeV. The peak in the invariant mass distribution occurs at E — |7| 2 /(2/i), where 
ItI = V^7rc + 7im- The value at the peak is proportional to (|7| + 7i m )~ 1 - The full width at 
half maximum is 2(2|^y | + 7im)[(|7l + 7im)(3|7| + 7im)] 1/2 7 '/x. 

VI. THE X LINE SHAPE 

The X is observed as a peak in the invariant mass distribution of its decay products, 
such as J/ipn + n^. Its mass and width are extracted from that invariant mass distribution. 
For instance, the Belle collaboration obtained their value for the mass and the upper bound 
on the width by fitting the J/ipTT + ir~ invariant mass distribution in B + — > J/ip 7r + 7r~ K + 
near the D°D*° threshold to a resolution-broadened Breit-Wigner function on top of a 
polynomial background. The shape of the invariant mass distribution of the decay products 
of the X is called the line shape. The resonant interactions in the D D*°/D*°D system can 
significantly modify the line shape, so it need not have the conventional Breit-Wigner form. 
In this section, we compute the line shape of the X in short- distance decays of the X. To 
be definite, we consider the production process B — > HK, where H is the hadronic system 
consisting of J/ipTC + n~ with invariant mass near the D°D*° threshold. However our results 
on the line shape will apply more generally to any short-distance production process for X 
and any short- distance decay mode of X. 

The separation between the long-distance scale l/\a\ and all the short-distance momen- 
tum scales in the decay B — > HK can be exploited through a factorization formula for the 
T-matrix element: 

T[B -> HK] = Abort [B -> (DD*)° + K\ x A(E) x A short [(DD*)° + -> H\. (54) 

There is an implied sum over the spin states of the D*. The long-distance factor A(E) 
depends on the complex- valued scattering length a and is given in Eq. (35). Its argument 
E is the difference between the invariant mass M of the hadronic system H and the D°D*° 
threshold, as given in Eq. (51). The short- distance factor Ahort associated with the initial 
state is the same one that appears in the factorization formulas for B — > XK in Eq. (42) 
and for B — > D°D*°K in Eq. (48). The short- distance factor A s hoit associated with the final 
state is the same one that appears in the factorization formulas for X — > H in Eq. (32) and 
for D°D*° -> H in Eq. (34). 

The factorization formula in Eq. (54) can be motivated diagrammatically by separating 
the loop integrals in the decay amplitude according to whether the virtual particles are off 
their energy shells by less than or by more than A 2 /(2/i). There are terms in the decay 
amplitude for B — > HK that are enhanced by a factor of am n when the invariant mass of 
H is near the D°D*° threshold. These terms can be represented by the Feynman diagrams 
in Fig. 7 and can be expressed in the form 

T[B^HK] = A {A) [B -> (DD*)lK]L(A,E)A(E)L(A,E)A {A) [(DD*)l -> H]. (55) 

The factors A^ A \ which are represented by dots in Fig. 7, are amplitudes in which all virtual 
particles are off their energy shells by more than A 2 /(2/i). The factors of L(A,E), which is 
given in Eq. (9), are the amplitudes for the propagation of the D and D* between contact 
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FIG. 7: Feynman diagrams for B — > J/ipir + ir K that are enhanced near the D°D*° threshold by 
a factor of am^. The blob represents the geometric series of diagrams given in Fig. 1. 



interactions. The condition \E\ <C A 2 /(2/i) implies that the \f—2[iE term in L(A, E) can 
be neglected. The resulting expression for the T-matrix element is the factorization formula 
in Eq. (54), with the short-distance factors A^ort given in Eqs. (44) and (38). 

The factorization formula for the T-matrix element in (54) implies a factorization formula 
for the invariant mass distribution for the hadronic system H near the D°D*° threshold. 
If the hadronic system consists of particles with momenta pi and invariant mass M, the 
factorization formula is 

<W : [B+^K + H] 



dM 

y \A(E)? x J |Ah„„[(M>*)° + - H\\ 2 (2kYS 1 (Ph - E,P.) U (2 ^g £< ■ (56) 

For M near the D°D*° threshold, the only significant variation with M is through the long- 
distance factor |^4.(i?)| 2 , where E is the energy defined in Eq. (51). If the complex scattering 
length is parameterized as in Eq. (12), the long-distance factor is 



/ mi 

,2 I. ,2 



E > 0. (57b) 



((2^)V2 + 7im )2 + 7 2 e 

This factor gives the line shape of the X. Note that for E > 0, the line shape does not 
depend on the sign of 7 re . However for E < 0, the line shape is completely different for 
7 re > and 7 re < 0. 

In the case 7 re > 0, the peak in the invariant mass distribution occurs below the D°D*° 
threshold by the amount 7 j? e / (2/x). The X line shape is illustrated in the upper panel of Fig. 8 
for 7im = 10 MeV and for three positive values of 7 re : 10, 20, and 30 MeV. If 7im < 7 re , the 
full width of the peak at half maximum is 2 7re7im / /i. The line shape for E < is symmetric 
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FIG. 8: The X line shape in a short-distance decay channel, such as J/tp ir + ir~ , for 7; m = 10 MeV 
and for various positive values (upper panel) and various negative values (lower panel) of 7 rc /7im- 
The horizontal axis is the difference E = M — (m^o + m D *o) between the invariant mass M and 
the D°D*° threshold. 



about the peak as a function of ji?! 1 / 2 but not as a function of E. If 7 im <C 7 re , the line 



shape in Eq. (57) is sharply peaked at E = - 
function: 

47T 3 7 r 



7^/ (2/x) and it can be approximated by a delta 



\A(E)\^^5(E + 1 l/(2f,)). 

H' lim 



(58) 



Note that the condition 7 im <C 7 re is equivalent to T x <C E x /k- 

In the case 7 re < 0, the peak in the invariant mass distribution occurs at the D°D*° 
threshold. The X line shape is illustrated in the lower panel of Fig. 8 for 7 im = 10 MeV 
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and for three negative values of 7 rc : —10, —20, and —30 MeV. The line shape has a cusp 
at E = 0. Bugg has proposed that the X can be identified with this cusp at the D°D*° 
threshold [22]. The normalization is the same in the upper and lower panels of Fig. 8. Note 
that the area under the cusp in the lower panel of Fig. 8 is much smaller than the area under 
the resonance in the upper panel for the same values of 7 im and |7 rc |. Thus a cusp seems 
less likely as an interpretation for the X(3872) than a resonance, although a quantitative 
analysis would be required to rule out that possibility. 

The integral over all energies of the line shape of a conventional Breit-Wigner resonance 
is convergent. In contrast, the integral of the line shape in Eq. (57) diverges logarithmically 
as the endpoints E m [ n and E m&x of the integral increase in magnitude. This follows from 
the fact that the line shape in Eq. (57) decreases as 1/\E\ for (2/j,\E\) 1 ^ 2 ^> |7 re |,7im- That 
expression for the line shape is of course only accurate for \E\ lower than A 2 /(2/i) ~ 10 
MeV, where A ~ m n is the natural momentum scale for low-energy DD* scattering. Thus 



E„ 



the logarithmic dependence on E min and E max holds only for \E mh 
convenient to define p min and p max by E min = -p 2 min /(2fi) and E max = +p max 
integral of the factor in Eq. (57) reduces in the limit p m m,Pm&x |7rc|,7im to 



< A7(2/x). 
2 •(->//)• 



/' 

J E, 



\A{E)\ 2 dE 
where f(x) is the function 



Att 2 



bPminPmax . ^"7re 
g-*— — + 



Tre Tii 



7ir 



#(7re) - /(7re/7 im J 



f(x) = x arctan(l/rr) + (l/x) arctan(x) 



It is 
The 



(59) 



(60) 



This function has a limited range, varying from 1 at x — and x = ±oo to ir/2 at x — ±1. 

The factorization formula for the invariant mass distribution of H in the case 7 rc > 
has important implications for measurements of the branching fractions of X. Since the two 
short- distance factors in Eq. (56) are insensitive to E, we can set M to m D o +m D *o or to m x 
in those factors. The short-distance factor associated with the decay of the B + reduces to 
F[B — > XK]/(2ti\Ax\ 2 )- If 7rc > 0, the short-distance factor associated with the formation 
of H reduces to T[X — > H]/\Ax\ 2 - Thus the differential decay rate in Eq. (56) reduces to 



dY 

dM 



[B -> HK] = T[B -> XK] Br[X -> H] 



Tx\A(E)f 
2ttUv| 4 



(61) 



If the product of T[B — > XK] and Br[X — > H] is measured by integrating dT /dM over the 
energy interval from E min to E max with | E min \ , E max ^> Ex,Fx, it will be in error by the 
factor 



2ir\A 



\A{E)\ 2 dE 



1 + 



a 

16E XJ 



x 



1 + log 



\E n 



1/2 F l/2 



E x + T 2 X /(16E X ) 



- f{AE x /V x ) 



x 



AuE 



x 



.(62) 



The error would cancel in the ratio of the branching fractions for any two short-distance 
decay modes of X. The error would not cancel in the ratio of the branching fractions for 
a short-distance decay mode of X and one of the long-distance decay modes D°D°tt and 
D°Z)°7. This effect should be taken into account in analyzing the decays of the X(3872). 
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VII. SUMMARY 



If the X(3872) is a loosely-bound S-wave molecule corresponding to a C = + superposi- 
tion of D°D*° and D*°D°, the scattering length a in the (DD*)° + channel is large compared 
to all other length scales of QCD. The decays of the X implies that the large scattering 
length has an imaginary part. It can be conveniently parameterized in terms of the real 
and imaginary parts of 1/a as in Eq. (12). The binding energy E x and the width T x are 
expressed in terms of those parameters in Eqs. (14) and (13b). 

The large scattering length can be exploited through factorization formulas for decay rates 
of X. For short- distance decay modes that do not proceed through the decay of a constituent 
D* of the X, the long-distance factor in the factorization formula is proportional to l/\a\ 
and is given in Eq. (41). If a partial width of the X is calculated using some model with a 
specific binding energy for the X, the factorization formulas can be used to extrapolate the 
prediction to other values of the binding energy and to take into account the width of the 
X. 

The large scattering length can also be exploited through factorization formulas for pro- 
duction rates of X, D°D*° near threshold, D*°D° near threshold, and decay products of X 
with invariant mass near the D°D*° threshold. The long-distance factor in the factoriza- 
tion formula for production rates of X is proportional to l/\a\. For production of D°D*° 
and D*°D° near threshold, the factorization formula implies that the dependence on the 
invariant mass is through the factor y/2 fiE\A(E)\ 2 , where E is the invariant mass relative 
to the D°D*° threshold and |^4(E)| 2 is given in Eq. (53). The peak in the invariant mass 
distribution is above the threshold by the amount Ex + T x /(16Ex)- The line shape of the 
X can be measured through the invariant mass distribution of its decay products. In the 
case of short-distance decay modes, the factorization formulas imply that near the D°D*° 
threshold, the shape of the invariant mass distribution is given by the factor |^4(E)| 2 in 
Eq. (57). If Re(a) < 0, the distribution has a cusp at E — 0, as shown in the lower panel 
of Fig. 8. If Re(a) > 0, the distribution has a peak at E — —Ex, as shown in the upper 
panel of Fig. 8. In contrast to a Breit-Wigner resonance, the integral over the line shape 
is logarithmically sensitive to the endpoints of the integration region. This effect should be 
taken into account in analyzing the production and decay of the X. 
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